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Abstract 

This manual describes version 2.2 of the Monte Carlo event generator UNIBAB 
for large angle Bhabha scattering at LEP and SLC. UNIBAB implements higher order 
electromagnetic radiative corrections and the effects of soft photon exponentiation 
in a photon shower approach. Weak corrections are included through the use of an 
electroweak library. 
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Program Summary: 

• Title of program: UNIBAB, Version 2.2 (December 1995) 

• Program obtainable from: 

Ftp : / /crunch . ikp . physik . th-darmstadt . de/ pub/ anlauf /unibab 

• Programming language used: FORTRAN-77 

• Computer/Operating System: Any with a FORTRAN-77 environment 

• Number of program lines in distributed program, including test data, 
etc.: ^ 9000 (Including comments) 

• Keywords: radiative corrections, large angle Bhabha scattering, multiphoton ra- 
diation 

• Nature of physical problem: Higher order leading logarithmic QED radiative 
corrections to large angle Bhabha scattering (e"'"e~ — > e"'"e~ + n'j) including weak 
corrections at Z energies and beyond 

• Method of solution: Monte Carlo event generation 

• Restrictions on the complexity of the problem: The program assumes domi- 
nance of the s-channel Z exchange contribution and does not include the interference 
between initial and final state radiation. Therefore, the program is primarily de- 
signed to be used in the vicinity of the Z peak and only for not too small and not 
too large scattering angles (10° < 9* < 170°). With somewhat limited accuracy, the 
program is also useful for LEP2 energies. For cross checks, the contributions of the 
t-channel diagrams may be switched off, so that the program effectively simulates 
e"'"e~ — > yU'''/i~. However, the latter mode is not very efficient. 

• Typical running time: Strongly dependent on the energy and the chosen cuts. 
The test run took approximately 55 CPU seconds on a DEC Alpha 3000-600 running 
Digital Unix, 65 CPU seconds on an SGI Challenge running IRIX, 80 CPU seconds 
on an HP 9000/735 running HP/UX, 355 CPU seconds on an IBM RS/6000-520 
running AIX, and 485 CPU seconds on an Intel 486DX50 running Linux. 
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1 Introduction 



The study of electroweak physics in e^e~-colhsions close to the Z resonance has proven 
to be extremely fruitful The precision measurements of the Z parameters are very 

powerful in constraining the as yet unknown electroweak parameters and allow to put 
stringent limits on potential physics beyond the standard model. 

In the context of these precision measurements, small angle Bhabha scattering plays 
a major role because all LEP/SLC experiments use it for their luminosity measurements. 
Consequently, from the beginning much effort has gone into the construction of semi- 
analytical programs 0] and Monte Carlo event generators for Bhabha scattering in 
the small angle regime. With improved statistics, systematic errors begin to dominate, to 
which the theoretical uncertainties in the calculation of cross sections give an important, 
by now even sometimes dominating, contribution. At the level of a required precision of 
the order of 0.1%, the inclusion of higher order electromagnetic corrections is indispensable 
and is taken care of in state-of-the-art calculations and Monte Carlo programs, as reviewed 
in [0,i. 

In contrast, the situation for large angle Bhabha scattering is less favorable. In the 
vicinity of the Z peak, s-channel Z exchange and t-channel photon exchange are of similar 
importance, thus the calculation of radiative corrections is more involved than in the case 
of fermion pair production or small angle Bhabha scattering. For typical experimental 
cuts, the radiative corrections to large angle Bhabha scattering are numerically much 
larger than to small angle Bhabha scattering, and higher order electromagnetic corrections 
are of much higher importance. Nevertheless, large angle Bhabha scattering at LEPl and 
SLC in principle provides a sensitivity to the parameters of the electroweak sector of the 
Standard Model comparable to the other charged lepton final states, t^t~. In 

the energy regime of LEP2, the contribution of Z exchange to the Bhabha cross section 
is much smaller than at LEPl/SLC, hence it can be considered essentially as a general 
QED test. 

The state-of-the-art for large angle Bhabha scattering is reviewed in ||^. While semi- 
analytical calculations [0-0 provide rehable predictions for a predefined set of cuts, it is 
nevertheless experimentally desirable to provide a comparable precision in a full Monte 
Carlo event generator, which will then allow to impose arbitrary cuts, and to simulate 
the response of the detector to e"'"e~ e"'"e~ events. 

The present paper describes the implementation and status of the Monte Carlo event 
generator UNIBAB 2.2 that provides a level of precision comparable to semi-analytical 
programs. It is an update of the first published version 2.0 [|T^, and has been compared 
with other available Monte Carlo event generators (BHAGENE3 [|l|] and BHWIDE ^) as 
well as available semi- analytical calculations ^ both for LEPl and LEP2 energies. 

The distinguishing feature of UNIBAB from fixed order Monte Carlo event generators 
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is that it handles the leading logarithmic contributions 




« 6% (at LEP/SLC energies) (1) 

to the electromagnetic radiative corrections to all orders in a Monte Carlo parton shower 
algorithm, which includes the very important exponentiation of the soft photon contribu- 
tions automatically, as well as the contributions from multiple emission of hard coUinear 
photons. In addition, UNIBAB contains the full kinematics from multiphoton emission, 
since the radiated photons are generated explicitly. 



The first version of the Monte Carlo generator, UNIBAB89 [jT6|,|T^], was a "QED dresser" 
in the sense of i.e. it contained only the QED corrections to an effective Born cross 
section. It also suffered from several serious limitations (see appendix P). In version 2.0 
of UNIBAB, most of these limitations were overcome using the experiences with the Monte 



Carlo event generators KRONOS [|I^ and KROWIG ||19[ for deep inelastic scattering at HERA. 
Weak corrections are included through the use of an electroweak library. 

This write-up is organized as follows: In section ^ we outline the physics underlying 
the algorithms implemented in UNIBAB. The actual implementation is described in section 
^ The parameters controlling the execution of UNIBAB are discussed in detail in section ^, 
and the FORTRAN-77 interface is presented in section ^. Section ^ contains our conclusions. 
Distribution rules, the revision history, a listing of all external symbols and an example 
can be found in the appendices. 



2 Bhabha Scattering at High Energies 



In the structure function formalism |]7|,|20|-p2| the factorized expression for the differential 
cross section for the process 



e+(p+)e-(p_) e+(p;)e-(p'_) [7(^1) • ■ ■] (2) 



reads 

dV 



dtj-dt-dw+dw- 



I dx+da;_dy+dy_r(x+, Q')r(x_, Q')D{y+, Q')D{y^, 


d{tj^,t-,u+,u-) dt+dt_dM+du_ 



where a is the Born level cross section of the hard process, r{xi,Q'^) are the structure 
functions for initial state radiation, D{yi,Q^) are the fragmentation functions used for 
final state radiation, and is the factorization scale. The invariants are defined as 
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follows 



s = {p++pj)^ (4) 
s = {x+p+ + X-P-f = x+X-S (5) 



2 



t± = [p±-p'±) (6) 



/ \ 2 



t± = [x^p^-^ =^t^ (7) 

V y±/ y± 



u± = (p^ - P±) 
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u± = (x^Pt - — ] = —u± (9) 

V y±J y± 

where the electron mass has been neglected. In the absence of photon radiation we have 
of course t+ = t_ and = 

Here we restrict ourselves to the leading log approximation (LLA) where the structure 
functions r(x, Q^) and D{x,Q^) are identical. The structure functions D satisfy the 
evolution equation EI 



X 

Peeiz) = (11) 

with initial condition 

D{x,ml) =6{l-x) (12) 

An explicitly regularized version of (p!0D which is used in the Monte Carlo implemen- 
tation is given by 



X 

r l-e 



a 
2^ 



L 



2.1 The Algorithms for Multiphoton Radiation 
2.1.1 Initial State Radiation 



From (^-|^) it is obvious that the implementation of (|^) in a Monte Carlo event generator 
amounts to solving ([T^) by iteration and taking into account the energy loss in the hard 
cross section. 
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The initial state branching algorithm, which is implemented as a photon shower algo- 
rithm, has already been described in ref. 



1^ 



It has been taken over essentially unchanged, 
the only modification is that we require conservation of four-momentum at every emission 
vertex and allow the electron momentum to take on off-shell values after initial state radi- 
ation. The off-shell momenta of the electron and positron after initial state branching are 
then used as the input momenta for the subgenerator of the hard Bhabha cross section. 
In the present version, the factorization scale used for initial state radiation is fixed to 
Q' = s. 



2.1.2 Final State Radiation 

In the vicinity of the Z peak and for not too small scattering angles, the cross section 
is dominated by s-channel Z exchange. Therefore we use the following approach. If we 
neglect interference between initial- and final-state radiation, the radiative corrections in 
the final state will then be well approximated by the radiative corrections to Z decay. 

The differential decay rate for one-photon emission Z e"'"e~7 is related to the decay 
Z —>■ e"'"e~ at lowest order |^ : 

dr(3) ^ dr(2) dr(2) 

Ci{Xi, X2) ■ + C2{Xi, X2) 



dxidx2dn 



1,2 



df^i 



(14) 



where 



Ci(xi,X2) 
C2(Xi,X2) 



a 



XT 



27r(l-Xi)(l-X2) 



a 



27r(l-Xi)(l-X2) 



Here dF*^^'^) denotes the differential decay width of Z into e"'"e~ (e"'"e~7), df2i^2 refers to 
the direction of the outgoing e"*" (e~) in the Z rest frame, and Xi and X2 are the energy 
fractions of e"^ and e~ for the radiative decay. For a more detailed description we refer 
the reader to ||25|| . 

It is easily verified that, in the limit of emission of soft photons or photons (almost) 
coUinear with the outgoing fermions, eq. reduces to 

1 dF(3) a;? + x2 

(15) 



F(2) dxidxo 



(2pV ■ k){2p'_ ■ k) 



which exhibits the leading soft and coUinear singularities that have to be resummed to all 
orders. 

We choose the phase space boundary such that the leading logarithms are correctly 
reproduced. 



1 - ^3(1 -rf) < xi,2 < 1 - X'i'T] , e < X3 < 1 



(16) 



with 
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Xi + X2 + X3 



= 2 



The parameter rj is therefore related to the factorization scale associated with final 
state radiation: 

i2 



In 



In 



1 



;i7) 



In the case of multiphoton final state radiation approximations are necessary. Our 
main concern is the proper infrared and coUinear limit of the effective radiation matrix 
element, plus a good approximation for the photon that carries the largest transverse 
momentum relative to the outgoing fermions, while we want to neglect the interference 
between successive emissions, and assume that the radiated photons are widely separated 
in phase space. One can easily see that iteration of the single photon emission algorithm 
satisfies these criteria and reproduces the leading logarithms, provided the following vir- 



tuality ordering conditions are met. Denote by Xi 



X2 , 1 



1, . . . , n, the phase space 



variables for the case of emission of n photons. The phase space variables are then re- 
quired to satisfy 



mm 



1 



X 



(1) 



X2 



> mm 



J-2 



> 



> mm 



X 



(n) 



J-2 



1^ 



which are the analogue of the ordering conditions (15) resp. (17) of ref. [|T8|, except 
that they are reversed. This, together with the phase space (p!6D, guarantees the proper 
normalization of the leading logarithmic terms. 



PS 



C,ix?,xl^)+C2{x^>,x)^>) 





[/ 


n! 1 




IJ 


a 


ij! 1 


I- 






ij! ^ 







V 



21ni-? 

e 2 



(19) 



The algorithm for final state multiphoton radiation is thus defined as follows, i) 
Determine the number of photons Uf emitted from the final state according to ([T9|) . ii) 
In the case of single-photon emission, use distribution (|14|), with phase space limits (p!^). 
iii) For emission of multiple photons, iterate the algorithm for single photon radiation, 
taking into account the phase space conditions 



2.2 The Effective Born Amplitude 

The calculation of the effective Born amplitude follows the BHM/WOH approach as 
described in 0,^. The amplitude for e+e~ — > e+e" including virtual corrections can be 
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cast into a form close to the lowest order amplitude: 



A(e+e- ^ e+e-) = A^^) + Af^ + A^^^ + Af^ + {h 



'OX) 



(20) 



where A'^'^^ denotes the dressed photon, A^"^^ the dressed Z exchange amplitudes, and {box) 
the terms from box diagrams, which can be essentially neglected around the Z peak. 
The dressed photon exchange amplitude is written as 



while the Z exchange amplitude is given by 



[(1 + F^v)l, - FTlM ® [(1 + ^v^^)7'^ - FTl'l^] 



(21) 



^(Z) 



-geri"'(s))7^-(ae + Fi^)7M75 
+ F^^ + Q,Il^^(s))7^ - (ae + )7'^75' 



(22) 



Here n'^^(s) = S'^^(s)/(s + S'^'^(s)). For more details on the notation we refer the reader 
to the literature. 



2.3 QED Virtual Corrections 

Besides the weak virtual corrections that are easily included in the effective Born am- 
plitude (|20|) , we have to include the QED virtual corrections (vertex corrections, box 
corrections). In addition, we have to properly include the contributions from soft-photon 
emission beyond leading logs. 

The method of inclusion of nonleading QED corrections in the structure function 



formalism is well known (see e.g. ref. [^). However, to preserve a fully factorized form 
of the cross section similar to (|^), which is mandated by our implementation of initial- 
and final-state radiation, we refrain from introducing /^-factors and proceed as follows. 

The effective Born amplitude is multiplied by (helicity-dependent) form-factors, as 
given explicitly e.g. in the appendix of ref. 0]. The photon- mass singularities in the form- 
factors are explicitly canceled against the soft-photon contribution. The infrared- singular 
parts are factored out and absorbed into the structure functions and fragmentation func- 
tions, thus uniquely determining an infrared-finite nonleading correction. 

The cross section determined by this matching procedure in principle fully reproduces 
the 0{a) corrections in the soft limit. However, since the actual photon shower algorithms 
do not yet include initial-final interference, we have to drop the infrared-finite pieces of 
the QED boxes and corresponding terms from the soft-photon corrections for a consistent 
treatment . 
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2.4 Implementation of the Bhabha Cross Section 



The differential cross section for tlie liard subprocess is implemented in a standard way, 
using importance sampling techniques. As a majorant function for the differential cross 
section, we use (r = —i/s = (1 — cos^cms)/2) 



-(.,r)=A(.) + -^ 



(23) 



The functions A{s) and -B(s) are chosen such that the peaking behaviour of the full cross 
section for small angles as well as the resonant behaviour for large angles in the vicinity 
of the Z peak is efficiently reproduced: 



where 



A{s) = 
B{s) = 

Xl{s) = 



A'Ka{sf 



Xlis) 



Mr 



+ 



M|r| 



Mil 



{s - Mir + MiVl 



(24) 
(25) 

(26) 



Here a{s) is the running QED coupling constant, and v^. and Oe denote the vector and axial 
vector couplings of the electron to the Z. The parameter C 0.1 has been empirically 
adjusted so that da/dr is a majorant of the effective Born cross section. 

The angular distribution is generated by a combined mapping and rejection algorithm 
with mapping function 



^1 

B{s) T 

and corresponding event weight. 



X(t) = X(rmax) ■ P + X(rmm) • (1 - p) , 



w 



do"(s, r) 



(27) 



(28) 



d(7(s,T) ' 

where da is calculated using the effective Born amplitude, as described in the previous 



sections. 



3 Implementation of UN I BAB 2 . 2 

Like almost all Monte Carlo event generators, UNIBAB is divided into three parts: initial- 
ization, generation, termination. These are described in this section. 
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3.1 Initialization of UNIBAB 



The initializations in UNIBAB are used for computing the value of variables that will 
be used frequently during event generation. The primary example is the calculation of 
electroweak couplings, widths, and masses from a set of input parameters. To this end, 
the generic initialization routine ubinit calls the subroutine ubigsw which contains a 
simple interface to the initialization routine of the electroweak library. 

Furthermore, internal steering parameters are derived from the cuts specified by the 
user. Finally, a standard /hepevt/ initialization record is written [|T|, which can be read 
by the analysis program. 



3.2 Event Generation 

The routine ubgen produces an event on every call. The four momenta are written to a 
standard /hepevt/ event record, where they can be read by user supplied analyzers 
(the default configuration uses hepawk ||2B|). See section ^ for details on UNIBAB's use of 
/hepevt/. 

After generation of the initial state radiation by the branching routine ubbini, a raw 
Bhabha event corresponding to the approximate cross section a (ubxtot) is generated at 
reduced lepton energies with the corresponding angular distribution (ubgt). This event 
will be given a weight according to the ratio of the full differential cross section including 
weak corrections (ubxdif ) to the approximate one (ubxtri). 

Next, the routine ubgppr constructs the four- vectors of the leptons after the hard 
scattering in the effective c.m.s. frame. Finally, after generation of the final state radiation 
by the branching routine ubbf in, the event will be accepted in ubgacc if it passes the 
given experimental cuts, otherwise it will be rejected and the above algorithm repeated. 



3.3 Termination 

The cross section for the generated events, and its error, are obtained from the standard 
formulae 

, , r / # of successful trials , , 

A(T,„,(s) = max{a{s')}- (30) 

s'>so 



(total # of trials — # of successful trials) ■ # of successful trials 
^ (total # of trials)^ 

(where sq is the minimal invariant mass squared of the outgoing electron positron pair as 
determined from the given cuts) and placed into /hepevt/ for inspection by the analysis 
program, which might use it to normalize its histograms at this point. 
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4 Parameters 



The parameters controlling UNIBAB version 2.2 are summarized in table |^. 
4.1 Elect roweak Parameters 

Starting with UNIBAB version 2.0, the electroweak corrections are included through calls 
to an electroweak library. 

The electroweak library of UNIBAB 2.2 is based on the library that is distributed with 
ALIBABA 2.0 [0| and assumes the minimal standard model. This library is based on 
ref. which uses the on-shell renormalization scheme [^]. The input parameters are 



the masses of the Z (masslz), the top-quark (masslt) and the Higgs (masslh). During 
initialization, the mass of the W (masslw) as well as Fz (gammlz) and sin^ 9w (sin2w) 
are calculated from these parameters. As of version 2.2, the electroweak library has been 
updated to include QCD-corrections to the Tz and the leading m^-corrections to the 
p-parameter (see e.g., |@]). 

UNIBAB has several switches that control the "physics" entering the hard scattering 
cross section. The first, tchann, may be used to switch off the contributions of the t- 
channel diagrams to the cross section, which is useful for the application of the so-called 
"t-channel subtraction" to experimental data, and to compare with programs that describe 
the process e+e~ —>■ fi'^fj^~ including higher order QED corrections, like e. g. KORALZ pO |. 



With tchann set to . false . , the program effectively simulates the process e"^e~ /^^/^~; 
however, this mode of operation is not very efficient. 

The switches qedvtx and qedbox control the inclusion of the finite pieces of the QED 
vertex and box corrections. The QED vertex corrections are enabled by default, however, 
as explained in section \2,.'S[ the box corrections are disabled. The box corrections may 
be switched on to estimate the size of the neglected initial-final interference effects, but 
should not be used for physics simulations, since initial-final interference is not included 
for hard photon emission. 

The next switch, weak, may be used to switch off the pure weak corrections in the 
electroweak library. According to standard terminology, this includes the propagator 
corrections to the exchanged photon and Z. It should be noted, however, that Tz as 
well as sin^ 9]y are computed from the input parameters even if the weak corrections are 
switched off. A final switch, boxes, controls the inclusion of the pure weak box diagrams. 
The relative contribution of these diagrams is of the order of at most 0.1-0.2% in the 
energy range of LEP, whereas their computation may be quite time consuming. 

4.2 Control of QED corrections 

The photon shower algorithms are controlled by the following switches: 
• isrtyp: Key for initial state radiation 
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Variable name 


semantics 


Default value 


ahpla 


1 / q;qed 


137.0359895 


massle 


/"c 


0.51099906 • lO"'^ GcV 


masslz 


Mz 


91.1887 GeV 


masslt 


nit 


174.0 GeV 


masslh 


Mh 


300.0 GeV 


alphas 


as(Mz) 


0.124 


ebeam 


e^ beam energy 


46 GeV 


epol 


e" beam polarization 





ppol 


e"*" beam polarization 





ctsmin 


minimum cos^* 


-0.9 


ctsmax 


maximum cos^^* 


+0.9 


acocut 


maximum e^e^ acoUinearity angle 


20° 


ecut 


minimum outgoing e^ energy 


20 GeV 


evisct 


minimum invariant mass of final state 





nevent 


Number of events 


1000 


tchann 


Switch for t-channel diagrams 


.true. 


qedvtx 


Switch for QED vertex corrections 


.true. 


qedbox 


Switch for QED box corrections 


.false. 


weak 


Switch for weak corrections 


.true. 


boxes 


Switch for weak box diagrams 


.true. 


isrtyp 


Key for initial state radiation 


1 


f srtyp 


Key for final state radiation 


3 


epsiln 


Internal infrared cutoff 


10-5 


rseed 


Random number seed 


54217137 


errmax 


maximum error count 


100 


verbos 


verbosity 





runid 


run identification 




stdin 


standard input 


5 


stdout 


standard output 


6 


stderr 


standard error 


6 



Table 1: Parameters controlling UNIBAB. These range from standard model parameters 
(Mz, rrit, Mh, . . .), over experimental cuts to purely technical parameters like input/output 
units. 
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Initial state radiation disabled. 

1 Initial state radiation enabled, factorization scale = s. 



• f srtyp: Key for final state radiation 

Final state radiation disabled. 

1 Final state radiation enabled, factorization scale Q'^ = s (c.f. ([T7|)). 

2 Same, but with Q'^ = —i ■ e, so that ln(l/?7) = ln(— t/me) 

3 Same, but with = si/u 

The different choices for the factorization scale for final state radiation, which are of course 
completely equivalent at the leading logarithmic level, can be used to estimate the effect 
of the missing initial-final interference effects in UNIBAB. 



The remaining QED corrections (vertex, boxes) are described in subsection ^A. 



The former switch bstyle (UNIBAB versions 2.0 and 2.1) is no longer avaible, it has 
been superseeded by i srtyp and f srtyp. 



4.3 Cuts 

The region of phase space where UNIBAB generates events is controlled by specifying 
kinematical cuts. For the sake of simplicity and efficiency, there are five parameters: 

• ctsmin: minimum cos 6'* (see below), 

• ctsmax: maximum cos 6**, 

• acocut: maximum acoUinearity angle of the outgoing leptons, 

• ecut: minimum energy of either of the outgoing leptons, 

• evisct: minimum invariant mass of final state. 



The meaning of the angle 6* is the following. Let 6+ and 6'_ be the angle betweeen the 
outgoing positron and electron and the incoming electron beam, respectively. The angle 
9* is then defined by 



cos 




-+71 




COS 




_ — TT 





If all emitted photons are emitted strictly coUinearly with respect to the radiating lepton, 
9* equals to the scattering angle in the partonic subsystem. 

Since UNIBAB is designed for the calculation of QED corrections in the LLA in the large 
angle regime, the angular cuts are restricted by the requirement that the factorization scale 
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be the c.m.s. energy s, so that terms of the order a/vr ■ log(— t/s) and a/vr ■ log(— u/s) 
remain small compared to the leading logarithmic terms o/tt ■ log[s/ml). We require (cf. 
eqs. TO) 

a/vr ■ log(-£/s)| , \a/n ■ log(-n/s)| < 1% . (32) 
Therefore, the angular cuts in UNIBAB are restricted to the range 

- 0.985 < cos < cos 9*^^^ < 0.985 , (33) 

which corresponds to 

10° < e* < 170° (34) 

While the parameter ecut has a physical meaning, the cut evisct has a technical 
interpretation and is to be understood in the following way. Due to the implementation, 
each radiated photon may be associated with either initial state or final state. Thus 
one can define an "invariant mass of the final state". It should be noted, however, that 
this quantity is unphysical and should be used only for optimization of the Monte Carlo 
speed for a given set of physical cuts. An example is the calorimetric measurement of the 
Bhabha cross section, where one does not discriminate between the outgoing lepton and 
an accompanying photon if the latter lies within a narrow cone of the former, but where 
one cuts on the sum of their energies. 

Although most of the above-mentioned cuts have a simple physical meaning, their im- 
plementation in the Monte Carlo generator is not without problems. For strictly collinear 
photon emission, the relations (Q-^ between laboratory variables and the variables for 
the hard scattering subprocess allow a straightforward derivation of the Monte Carlo's 
internal steering parameters from these cuts in a strict manner. If one allows for finite 
angles between the photon and the radiating lepton, these relations get modified, and 
one has to weaken the internal cuts; in fact, one cannot even derive these parameters any 
more in a strict manner, since the multi-particle phase space becomes too complicated. 
In the actual implementation, we have used empirical relations to derive these internal 
steering parameters from those valid for collinear emission, so that the contribution to 
the systematic error on inclusive observables like cross sections is well below the overall 
precision, which is of the order of 1%. For less inclusive observables, we would recommend 
that the user convinces himself by a variation of the cuts that the generated distribution 
does not change in the region of interest. 

It should be noted that setting the lepton energy cut ecut to very small values or 
setting the acoUinearity cut acocut close to 180° makes the Monte Carlo generator very 
slow. The latter also implies that non-logarithmic terms, which are not implemented in 
UNIBAB, become more important. 
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* unibabappl . f 



call ubdcmd 



Cinit') 



initialize the generator 



call ubdcmd 



('generate 10000 



generate 10000 events 



call ubdcmd 



Ccl 



ose ' ) 



cleanup 



Figure 1: FORTRAN-77 interface 



4.4 Monte Carlo Parameters 

The remaining, more technical Monte Carlo parameters should be almost self explaining. 
Since our photon shower algorithm automatically includes soft photon exponentiation, 
the results will not depend on the value of the internal infrared cutoff epsiln (which is 
measured in units of the beam energy), provided it is kept well below the experimental 
energy resolution. However, it is not advisable to set it many orders of magnitude lower 
than the default value, because this may result in too high photon multiplicities which 
will overflow internal tables. 



UNIBAB version 2.2 provides two application program interfaces on different levels. The 
higher (preferred) level consists of the command interpreter ubdcmd which accepts com- 
mands in form of character* (*) strings. This driver communicates with the analyzer 
hepawk |2^ by default. The lower level consists of two FORTRAN-77 subroutine calls: 
ubpsrv and unibab. 

5.1 Higher Level Interface 

The simple commands understood by ubdcmd are (here keywords are typeset in typewriter 
font and variables in italics; vertical bars denote alternatives) 

• initialize 

Force initialization of UNIBAB and write an initialization record into the /hepevt/ 
event record, which should trigger the necessary initializations in the analyzer. 

• generate [n] 



5 FORTRAN-77 Interface 
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Generate nevent events and call hepawk to analyze them. If the optional parameter 
n is supplied, nevent is set to its value. 

close 

Write a termination record to /hepevt/, which should trigger the necessary cleanups 
in the analyzer. 

statistics 

Print performance statistics (this is usually only useful for the UNIBAB developers, 
who are tuning internal parameters). 

quit 

Terminate UNIBAB without writing a termination record. 
exit|bye 

Write a termination record and terminate UNIBAB. 
set variable ival\rval 

Set physical or internal parameters. See the table |l] for a comprehensive listing of 
all variables. For example, the command 'set ahpla 128.0' will set the QED fine 
structure constant to 1/128. 

print variahle\all 

Print the value of physical or internal variables. Specifying the special variable all 
causes a listing of all variables known to UNIBAB. 

debug|nodebug flag 
Toggle debugging flags. 

testran 

Test the portability of the random number generator. We use a generator of the 
Marsaglia-Zaman variety |31| , which should give identical results on almost all ma- 
chines. 

banner 

Print a string identifying this version of UNIBAB. 
echo message 
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* ubdcmd.f 

subroutine ubdcmd (cmdlin) 
character* (*) cmdlin 

else if (cmdlin. eq. ' init ' ) 
call unibab (0) 
call hepawk ('scan') 
else if (cmdlin. eq. 'generate' ) 
do 10 n = 1, nevent 
call unibab (1) 
call hepawk ('scan') 
10 continue 

else if (cmdlin. eq. 'close') 
call unibab (2) 
call hepawk ( ' scan ' ) 

else 
end 

Figure 2: Event generation loop 

Print message on standard output. 

Unique abbreviations of the keywords are accepted, i.e. 'g 1000' generates 1000 
events. The tokens arc separated by blanks. Blank lines and lines starting with a # 
are ignored and may be used for comments. For portability, only the first 72 characters 
of each line are considered. 

On UNIX systems UNIBAB reads default startup files .unibab in the user's home 
directory and the current directory, if they exist. 

5.2 Lower Level Interface 

The subroutine unibab (code) has a single integer parameter. The parameter code is 
interpreted as follows: 

• 0: initialize the generator and write an initialization record to /hepevt/. 

• 1: generate an event and store it in /hepevt/. If UNIBAB has not been initialized yet, 
the necessary initializations are performed, but no initialization record is written. 

• 2: perform final calculations and write the results to /hepevt/. 



initialize UNIBAB 
print initialization record 

generate an event 
analyze the event 

get total cross section, etc. 

finalize analysis 
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Figure |^ displays excerpts from a simplified version of ubdcmd which make the correspon- 
dences between the two levels of the FORTRAN-77 interface explicit. 

UNIBAB's parameters can be accessed on the lower level by the subroutine ubpsrv 
(result, action, name, type, ival, rval, dval, Ival). The parameter is speci- 
fied by its (lowercase) name in the character* (*) string name. The string action is 
either 'read' or 'write' corresponding to whether the parameter is to be inspected or 
modified. The type of the parameter ('int', 'real', 'dble', or 'Igcl') is returned 
in type, if action is set to 'read'; if action= 'write ' , the type must be specified in 
type. Depending on this type the value is passed in ival, rval, dval, or Ival, respec- 
tively. The following error codes will be returned in the string result: ' ': no error, 
'enoarg': invalid action, 'enoent': no such parameter, 'enoperm': permission denied, 
and 'enotype': invalid type. 

The protection scheme implemented with this parameter handling has been described 
TB| . Its main purpose is to guarantee consistency of user defined and computed pa- 
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rameters in the generation phase of the Monte Carlo. 
5.3 Additional information in /hepevt/ 

Because UNIBAB uses since version 2.2 the standard /hepevt/ event record internally, not 
only stable particles with isthep(i) = 1 will be present. Adapting the conventions of 



the HERWIG Monte Carlo we use the following status codes 

• 101: e~ beam, 

• 102: e"*" beam, 

• 103: e"'"e~ center of mass system, 

• 110: e+e~ hard scattering center of mass system, 

• 111: e~ before hard scattering, 

• 112: e"*" before hard scattering, 

• 113: e~ after hard scattering, 

• 114: e"*" after hard scattering. 

However, these entries have no physical significance and should never be used in any anal- 
ysis. An exception to this rule are the beam particles 101 and 102, which are convenient 
for defining the reference frame (positive z axis pointing in the direction of the incoming 
electron) and are used e.g. by hepawk for this purpose. Only the particles with status 
code 1 belong to the final state as predicted by UNIBAB. 
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6 Conclusions 



We have presented the new version 2.2 of the Monte Carlo event generator UNIBAB for 
Bhabha scattering at LEP/SLC. The distinguishing feature of UNIBAB is the inclusion 
of higher order electromagnetic corrections, including exponentiated soft photons, in a 
photon shower approach. In contrast to fixed order calculations which have to be expo- 
nentiated by hand, UNIBAB handles the multiphoton effects explicitly. 
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A Availability 

The latest release of UNIBAB is available from the authors upon request. In addition, 

UNIBAB is available by anonymous ftp from the directory 

ftp : //crunch . ikp . physik . th-darmstadt . de/pub/ anlauf /unibab 



It is nevertheless recommended to notify the authors, in order to be informed of future 
bug fixes and enhancements. 



UNIBAB is distributed in PATCHY/CMZ card format |3^, 35], but machine specific 
and plain FORTRAN-77 versions are made available on request. A more modern (auto- 
configuring) and self-contained version of the Monte Carlo generator will be made available 
in a future release. 



B Revision History 

Version 2.2, December 1995 

• Update of the electroweak library. 

• Improved final state photon shower algorithm. 

• Different choices of factorization scale for final state radiation. 

• QED virtual corrections. 

• Speed-up of electroweak corrections. 

• Fixed minor bug in initial-state photon shower. 
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Version 2.1, September 1994 

• Beam polarization. 

Version 2.0, June 1993 

• weak corrections included. 

• problems with energy non-conservation of old final-state radiation algorithm fixed. 

• more user-friendly: physical cuts replace "bare" Monte Carlo steering parameters. 

Version 1.3, September 1992 



• improved form of initial state radiation, as implemented in . 

• final state radiation. 

Version 1.2, July 1992 

Complete internal rewrite 

• /hepevt/ support (also used internally now), 

• double precision arithmetic (which is the native floating point mode on most modern 
machines) , 

• new driver program (using KRONOS [|18[ routines). 

The physics algorithms have not changed from version 1.1. The purpose of this release is 
to provide a usable baseline for further improvements. 

Version 1.1, September 1989 

Maintenance release. 

Version 1.0, July 1989 

First official release. 
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C Common Blocks and Subroutines 



To avoid possible name clashes with other packages, all external symbols exported by 
UNIBAB proper begin with the two letters UB, except for the routine unibab itself and 
the /hepevt/ common block. The names of the subroutines and common blocks of the 
provided electroweak library have been taken over unchanged. 

• Common Blocks: 

— /ubpcom/: Main parameter common block, holds all physical parameters. Ap- 
plication programs should access this common block through the ubpsrv rou- 



tine (cf. section p^A I- 

/ubcbrn/: holds the maximum Born cross section in the available energy in- 
terval. 

/ubcevt/ 

/ubcmsc/ 

/ubcsta/ 

/ubctri/ 



passes the particle momenta between subroutines, 
internal parameters. 

holds statistics on UNIBAB's performance, 
lookup table for driver keywords. 



Primary entry point: 



— unibab: main entry point, described in section |5l2 . 
Branching: 

— ubbini: initial state branching for one lepton. 

— ubbf in: final state radiation for one lepton. 

Cross sections: 

— ubxdif: differential cross section d(T/dcos6' (in pbarn). This function contains 
a very simple interface to the electroweak library. 

— ubxtri: approximate differential cross section used for mapping. 

— ubxtot: integral of the approximate cross section used for mapping. 

— ubxcof : coefficients of the approximate cross section. 

— ubf f ac: QED vertex form factors and box corrections. 

Initialization: 

— ubinit: initiahze UNIBAB. 
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— ubigsw: setup standard model parameters. This subroutine also calls the ini- 
tialization routines of the electroweak library and transfers parameters between 
these two modules. 

— ubibmx: calculate the maximal Born cross section in the considered energy 
interval. 

— ubibn: negative Born cross section. 

— ubeeni: enters an initialization record into /hepevt/. 

Event generation: 

— ubgen: main event generation routine. 

— ubgppr: construct four momenta of hard subprocess in the effective cm. sys- 
tem. 

— ubgt: generate momentum transfer in the hard subprocess. 

— ubgacc: checks if current raw event satisfies physical cuts. 

— ubeent: enter a particle in /hepevt/. 

— ubenew: start a new /hepevt/ record. 

— ubenul: zero a /hepevt/ entry. 

Termination: 



ubclos 



ubeens 
ubstat 



write a final record to /hepevt/. 



enters a termination record into /hepevt/. 
print statistics on UNIBAB's performance. 

Driver program: 

— ubdriv: main driver program, calling the command loop. 

command loop, reading commands from the terminal or from files, 
execute a single command (cf. section |5.1| ). 
get the next string from the command line, 
get the next floating point number from the command line, 
get the next integer from the command line. 

Parameter management: 



ubdloo 
ubdcmd 
ubdlxs 
ubdlxd 
ubdlxi 



ubpsrv: inspect or modify parameters (cf. section |5.2|) . 

block data ubpini: initialize parameters with their default values. 
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ubrgen 
ubrtst 
ububoo 
ubulwr 
ubumin 
ubumsg 
ubuort 
ubupro 
ubut im: 



— ubpprn: print parameters. 

General utility routines: 

random number generator of the Marsaglia-Zaman type |HT 
selftest routine for ubrgen. 
boost a four-vector, 
convert string to lower case, 
minimize function, 
print a message. 



construct three orthogonal three-vectors for a given z direction, 
project a pair of four-vectors onto mass shell, 
check remaining CPU time. 

Keyword search (using the dynamic tries described in ||36|| ): 

— ubtins: insert a new keyword. 

— ubtlup: look up a (possibly abbreviated) keyword. 

— ubtnew: insert a new node into the trie. 

— ubtlen: calculate length of keyword. 

— ubtc2a: convert keyword from character* (*) to integer (*). 
Electroweak library entry points: 

— ewinit: initialize electroweak library. 

— eeeew: effective cross section including virtual corrections. 
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Test Run 



UNIBAB version 2.2 is distributed together with a sample command file and hepawk script, 
which are given below. To run this example, the user will need to link UNIBAB with the 
CERN library, because histogramming is done by HBOOK . 



The file sample .unibab is read from standard input (unit stdin, which is initialized 
to 5), and sample .hepawk is read from the file SCRIPT (i.e. under MVS from the file which 
has been allocated to the DDNAME SCRIPT and under UNIX from the file script or from 
the value of the environment variable SCRIPT). 

After UNIBAB has been successfully built and linked with hepawk on a supported UNIX 
system, running either maike check or make test should validate the build. 

sample .unibab 

Here is a simple UNIBAB command file, setting up parameters and generating 10000 events. 



# sample .unibab — sample UNIBAB command file 



# kinematical cuts (e+e-) 

set ctsmin -0.8 

set ctsmax 0.8 

set ecut 20 

set acocut 30 

set ebeam 48 



# run 
init 

gen 10000 

close 

quit 
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sample .hepawk 



This is a small hepawk script that plots the acoUinearity distribution for Bhabha events, 
for which the acoUinearity is larger than 0.01 radians. 

# sample .hepawk — sample HEPAWK analyzer for UNIBAB. 

BEGIN 
{ 

printf ("\nWelcome to the UNIBAB test:\n"); 

printf ( " ***************************\n\ii" ) ; 

printf ("Monte Carlo Version: '/>s\n" , REV); 

printf (" Run: °/.d\n" , RUN); 

printf (" Date: y.s\n\n" , DATE); 

lower_angle_cut = 20/DEG; 
upper_angle_cut = PI - lower_angle_cut ; 
minimum_acollinearity_cut = 0.01; 

h_acollinearity 

= bookl (0, "electron acoUinearity", 50, 0, PI/6); 

incut =0; # initialize counter 

} 
{ 

# Collect the outgoing electron and positron 
$electron = $positron = $NULL; 

for ((§1 in LEPTONS) 

if (@l:id == _pdg_electron) 

$electron = @l:p; 
else if (01: id == - _pdg_electron) 

$positron = @l:p; 

} 

lower_angle_cut <= angle ($electron, (§Bl:p) <= upper_angle_cut 
&& lower_angle_cut <= angle ($positron, @Bl:p) <= upper_angle_cut 
{ 

# sample the acoUinearity, but cut out Born-like events 
acoUinearity = PI - angle ($electron, $positron) ; 

if (acoUinearity > minimum_acollinearity_cut) 
{ 

incut ++; 

fill (h_acollinearity , acoUinearity) ; 
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} 

} 



END 
{ 



# Dump some numbers 
printf ("\nRESULTS:\n") ; 
printf ("********\n\n") ; 

printf ("generated events: generated cross section: 'hg nb\n", 

NEVENT, XSECT * le6) ; 
printf ("events within cuts: '/d, cross section: °/og nb\n" , 

incut, (incut/NEVENT) * XSECT * le6) ; 



printf ("\nHISTOGRAMS:\n") ; 
printf ( " ***********\n\n" ) ; 
scale (XSECT/NEVENT*le6) ; 
plot ; # plot the histograms 

printf ( " \ndone . \n" ) ; 

} 
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sample . output 

The following output should result from the input files above, modulo small roundoff 
errors. 

ubdcmd: message: Starting UNIBAB, Version 2.01/08, (build 951206/1241) 
hepawk: message: starting HEPAWK, Version 1.6 

Welcome to the UNIBAB test: 

Monte Carlo Version: v02.01 (Dez 06 00:00:00 1995) 
Run: 1035996352 
Date: Dez 06 13:07:00 1995 



RESULTS : 

generated events: 10000, generated cross section: .2484 nb 

events within cuts: 5904, cross section: .1467 nb 

HISTOGRAMS : 

lelectron acollinearity 
HBOOK ID = 1 DATE 06/12/95 
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lOIO 

OIOIOO I 

lOIOOO 01 
I 100 000000000000 lOOOIIO 
II IIIIOOII 001 



CHANNELS 10 
1 



1 2 3 4 5 

12345678901234567890123456789012345678901234567890 



CONTENTS 1. 1 11111 

*10** 2 02888001207543221111110101000000000000000000000000 

863040 19453955855746329290866974455566444423445334 

97779369787190039861961166622697994727229927292474 

40246609138263740965179798911952767111227743272837 

LOW-EDGE 1. 1111111111222222222333333333344444444455 

*10** 1 01234567890123456789013456789012345678901235678901 

00011223344556677888990011223344556667788990011223 

04948383727161605059493838272616150594948382727161 

07429630852964185207418630742963085296418529741853 



* ENTRIES 
.00 

* BIN WID 
.10 



11808 
,1047E-01 



* ALL CHANNELS = . 1467E+00 



* MEAN VALUE 



,1140E+00 



* UNDERFLOW = 



* R . M . S = 



done . 

ubdriv: message: bye. 
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